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We utilize nonequilibrium covariant transport theory to determine the region of vahdity of causal 
QQ ] Israel-Stewart dissipative hydrodynamics (IS) and Navier-Stokes theory (NS) for relativistic heavy 

. ion physics applications. A massless ideal gas with 2—^2 interactions is considered in a O+ID 

• Bjorken scenario, appropriate for the early longitudinal expansion stage of the collision. In the scale 

' invariant case of a constant shear viscosity to entropy density ratio rj/s fa const, we find that Israel- 

Stewart theory is 10% accurate in calculating dissipative effects if initially the expansion timescale 
exceeds half the transport mean free path To/Xtr,o ^ 2. The same accuracy with Navier-Stokes 
requires three times larger ro/Atr,o ^ 6. For dynamics driven by a constant cross section, on the 
other hand, about 50% larger To/Xtr,o <^ 3 (IS) and 9 (NS) are needed. For typical applications at 
' RHIC energies ^/SNN ~ 100 — 200 GeV, these limits imply that even the Israel-Stewart approach 

becomes marginal when ri/s > 0.15. In addition, we find that the 'naive' approximation to Israel- 
Stewart theory, which neglects products of gradients and dissipative quantities, has an even smaller 
range of applicability than Navier-Stokes. We also obtain analytic Israel- Stewart and Navier-Stokes 
solutions in 0-l-lD, and present further tests for numerical dissipative hydrodynamics codes in 1-1-1, 
^ ■ 2+1, and 3-l-lD based on generalized conservation laws. 

I. INTRODUCTION 



< 



J> ' The realization that shear viscosity is Ukely nonzero in general[l-3], and therefore the perfect (Euler) fluid 
CO i paradigm[4-7] of nuclear collisions at the Relativistic Heavy Ion Collider (RHIC) could have significant viscous 
' corrections [8], has fuelled great interest in studying dissipative hydrodynamics [9- 18]. Causality and stability 
problems[19] exhibited by standard first-order relativistic Navier-Stokes hydrodynamics[20, 21] steered most effort 
toward application of the second-order Israel-Stewart (IS) approach[22, 23]. 
OO However, unlike Navier-Stokes that comes from a rigorous expansion[24] in small gradients near equilibrium, the 

. IS formulation is not a controlled expansion in some small parameter (see Section II). Moreover, though causality 
is restored in a region of hydrodynamic parameters, the stability of IS solutions is not necessarily guaranteed [2 5]. 
Therefore it is imperative to test the applicability of the IS approach against a stable, nonequilibrium theory. 

In this work we perform such a test utilizing the fully stable and causal covariant transport approach[26-29]. 
We focus on the special case of 2 ^ 2 transport and a longitudinally boost invariant systcm[30] with transverse 
translational symmetry, i.e, 0-1-1 dimensions. Follow-up studies in higher dimensions, such as our earlier comparison 
d between transport and ideal hydrodynamics in 2-|-lD[8], will be pursued in the future. 

A similar study by Zhang and Gyulassy[27] compared kinetic theory and Navier-Stokes results. Here we compare to 
the causal IS solutions. In addition, we provide a series of tests and semi-analytic approximations that demonstrate 
the general behavior of IS solutions, which can be utilized to verify the accuracy of numerical IS solutions. 

The paper is structured as follows. We start with reviewing the relationship between hydrodynamics and covariant 
transport (Sec. II), then proceed to discuss the Israel-Stewart equations (Sec. III). The basic observables studied 
here are introduced in Sec. IV, while the main results from the hydro-transport comparison are presented in Sec. V, 
together with implications for heavy-ion collisions. Many details are deferred to Appendices A-D. We highlight 
here the generalized conservation laws derived in App. B, and the detailed study of Israel-Stewart and Navier-Stokes 
solutions in App. C utilizing numerical and analytic methods. 
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II. HYDRODYNAMICS AND COVARIANT TRANSPORT 



Hydrodynamics describes a system in terms of a few local, macroscopic variables [20], such as energy density s{x), 
pressurep(x), charge density n{x) and flow velocity u^{x). The equations of motion are energy-momentum and charge 
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conservation 

d^T^''{x)=0 , d^m{x) = o, (1) 
and the equation of state p(e,n). Ideal (Euler) hydrodynamics assumes local equilibrium in which case 

7t^,id = '^^«5(£,P,P,p) , N^R,id = i^^()) [<ii = (l,0)] (2) 

in the fluid rest frame LR. Extension of the theory with additive corrections linear in flow and temperature gradients [20] 

ST^I.^s = rjsiVu" + Vu" - Ia'^'O^Uc,) + (A'^'d^u^ , (3) 
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SN^s 



leads via (1) to the relativistic Navier-Stokes (NS) equations. (We use the Landau frame convention u^JT''^ = 
throughout this paper, i.e., the flow velocity is chosen such that momentum flow vanishes in the LR frame.) Here 
r]s{e,n) and C{e,n) are the shear and bulk viscosities, while Kq{e,n) is the heat conductivity of the matter. The most 
notable feature of NS theory relative to the ideal case is dissipation, i.e., entropy production. For consistency, the 
dissipative corrections (3)- (4) must be small, otherwise nonlinear terms and higher gradients should also be considered. 

It is crucial that the above hydrodynamic equations can indeed be obtained from a general noncquilibrium theory, 
namely on-shell covariant transport [21, 27-29]. For a one-component system the covariant transport equation reads 

p^'d^f{x,p) = S{x,p) + C[f,f]{x,p) (5) 

where the source term S specifics the initial conditions and C is the collision term. Throughout this paper we consider 
the Boltzmann limit [43] with binary 2^2 rates 

C[f, 9]{X, Pi) = jJJ {fs94 - fm) W"l2^34 S'^{pi+P2-P3-P4) (6) 
234 

where fi = f(x, pi) and = J d'^Pi/ (2Ei). For dilute systems, / is the phasespace distribution of quasi-particles, while 
the transition probability W — [1 / n) s{s ^ 4m'^)da / dt is given by the scattering matrix clcmcnt[21]. Our interest here, 
on the other hand, is the theory near its hydrodynamic limit, T4^ — > oo. In this case, "particles" and "interactions" do 
not necessarily have to be physical but could simply be mathematical constructs adjusted to reproduce the transport 
properties of the system near equilibrium[31]. The main advantage of transport theory is its ability to dynamically 
interpolate between the dilute and opaque limits. 

The Euler and Navier-Stokes hydrodynamic equations follow from a rigorous expansion of (5) in small gradients 
near local equilibrium 

fix, p) - fegix, p)[l + <Pix, p)] 101 « 1 , \p^'^^(|>\ « Ipf^dM/fe, , (7) 

and substitution of moments of the solutions 

N'^{x)= f^p^ f{x,p) , T^^-^{x)= f^p^p'^ f{x,p) , (8) 
J Po J Po 

into (1). The 0-th order </> = reproduces ideal hydrodynamics. The first order result is the solution to the linear 
integral equation 

P^dMx, P) = 2C[/eg, feg<l>Ns]{x, p) (9) 

and leads to the Navier-Stokes equations. 

Unfortunately the relativistic Navier-Stokes equations are parabolic and therefore acausal. A solution proposed by 
Mueller[32] and later extended by Israel and Stewart [22, 23] converts the NS equations into relaxation equations for 
the shear stress n^", bulk pressure 11, and heat flow q^. The dissipative corrections 

ST'^'^ = n'"' -UAf"" , SN^' = --^q^' (u^g^ = 0, u^tt^'^ = m^tt'^^ = 0) (10) 

dynamically relax on microscopic time scales r^(e,n), Tn(e,n), Tg{e,n) towards values dictated by gradients in flow 
and temperature. Causality is satisfied in a region of parameter space, however, stability is not guaranteed [25]. 
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More importantly, unlike the Euler and NS equations, the Israel-Stewart approach is not a controlled approximation 
to the transport theory (5). Instead of an expansion in some small parameter, it corresponds to a quadratic ansatz 
[23, 33] for the deviation from local equilibrium 



(11) 



where and are determined by the local dissipative corrections tt^'^, 11, and g^[44]. In contrast, the Chapman- 
Enskog solution (9) contains all orders in momentum. An evident pathology of the quadratic form (11) is that, in 
general, (j)G is not bounded from below and thus the phase space density becomes negative at large momenta (cf. (7) 
and (62)). Furthermore, the two approaches give different results not only for the relaxation times[21, 23], e.g., 

IS _ 

2p ' 



NS 



= 



(12) 



but also for the transport coefficients themselves. For an energy-independent isotropic cross section and ultrarela- 
tivistic particles (T ^ m) the difference is small[21], e.g.. 



NS 



AT 



0.8436— , r?f = — , (13) 

atr oatr 

but can be large for more complicated interactions. Here atr = J dilcm sin^ Ocmdcr /dClcm is the transport cross section 
(for isotropic, ctj^ = 2(Txot/S)- 

In the following Sections we analyze IS hydrodynamic solutions analytically and numerically, and test the accuracy 
of the IS approximation via comparison to solutions from full 2 — > 2 transport theory. 



III. ISRAEL-STEWART HYDRODYNAMICS AND BOOST INVARIANCE 



A. Israel-Stewart equations 



There seems to be some confusion regarding Israel-Stewart theory[22, 23] in the recent literature, therefore we start 
with reviewing the key ingredients. The starting point of Israel and Stewart (IS) is an entropy current that includes 
terms up to quadratic order in dissipative quantities [45] 



5" = !*^ 
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fin 



P 



T 



(14) 



(we follow the Landau frame convention). Here /x is the chemical potential, and Seq is the entropy density in local 
equilibrium. The coefficients {aj(e,n)} and {/3j(e,n)} encode additional matter properties that complement the 
equation of state and the transport coefficients. Most importantly, {Pi} control the relaxation times for dissipative 
quantities: 



(15) 



The entropy current and relaxation times in Navier- Stokes theory are recovered when all the coefficients are set to 
zero /3o = /3i = /32 = = ao = CKi (but as discussed previously, the IS and NS transport coefficients differ in general). 
The requirement of entropy non-decrease {d/^S^ > 0), which IS satisfy via a positive semi-definite [46] quadratic 

ansatz 



H2 

Td,,S^ = -r- 



KqT 



+ 



C KqT 2r]s 
leads to the identification of the dissipative currents: 

1„^., fffou"' 



> 



n 



c 



T 



- (3oDU + aadf.q'' - a'^q^Du^ 



-KqTA'"' 



Tn 



(|;) + \(l'^Tdx (^^\ + PiDq^ + aoV^H - aid^xu - aoUDu^ + oitta^ W 



a,- = 



2r?, 
d{l/T) 



1 „.™. f P2U^ 



(16) 

(17) 
(18) 
(19) 
(20) 



fj./T=const 
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Here D = u^dfj, and the () brackets denote traceless symmetrization and projection orthogonal to the flow 

A^^^-) = 1a''«A"^^^(A„^ + Ap^) - Ao.pA'-^ . (21) 

The new matter coefficients {aj(e, n)} are needed to describe how contributions from the q^^H and quir^" terms in (14) 
are split between the bulk pressure and heat flow, and the heat flow and shear stress evolution equations, respectively 
(in other words, a whole class of equations of motion generates the same amount of entropy - see Appendix A). 

Notice that the time-derivatives of heat flow, g'', and shear stress tensor, tt'^" are not expressed explicitly in (18)- 
(19) - instead, orthogonal projections to the flow velocity vector appear (cf. Eqs. (8a)-(8c) in [22]). Reordering 
the equations explicitly for the time derivatives gives rise to terms —u^qi,Du^ and —{tt^^u'^ + n^'^u^)Du\. There 
is therefore no need for a kinetic theory treatment [34] to derive these terms. They were missed in Ref. [12], but 
they are already present in standard IS theory as a trivial consequence of the product rule of differentiation and the 
orthogonality of the flow velocity and shear strcss/hcat flow. 

As we saw above, the Israel-Stewart procedure only determines the equations of motion up to nonequilibrium terms 
that do not contribute to entropy production. In kinetic theory, further such terms arise [23] when the vorticity 

= ^-A>^-A-P{dpu„ - d^un) (22) 
is nonzero. Including the vorticity terms, the complete set of evolution equations for the dissipative currents are: 

£»n = -^(n + cv^u^) (23) 



Bq^" = 



Po Po 

1 



T, 



- u^q^Du" (24) 



q L 



Pi Pi Pi Pi 

D-K^"" = -— fvr'"' - 27?V<''«''>) - (tt^^u^ + -k^''u^)Dux (25) 

-^TT"^ (vAt.^ + iPln|)-27r,<''c.'^)^ 

P2 P2 

We will refer to these equations as "complete IS" . If we ignore their tensorial structure, the equations have the general 
form 

X = - — {X-Xns) + XYx + Zx (26) 

TX 

for each dissipative quantity X, where Xj^s is the value of X in Navier-Stokes theory and Yy, Zx are given by 

the ideal hydrodynamic fields and dissipative quantities other than X. Therefore, Israel-Steward theory describes 
relaxation towards Navier-Stokes on a characteristic time tx, provided IYxIta: 1 and \Zx\tx jATivsj. 

In the last step of their derivation, Israel and Stewart neglect the first terms of the second lines, the ones with a 
factor 1/2 (this gives the Landau frame equivalent to (7.1a)-(7.1c) in [23]), because they expect to study astrophysical 
systems with small gradients {d'^u" + d''u'^\/T < 1, \d'^e\/{Te) < 1, ia^nl/(Tn) < 1, near a global (possibly rotating) 
equilibrium state. The neglected terms arc then products of small gradients and the dissipative quantities. We will 
refer to this approximation as "naive IS" [47]. In heavy ion physics applications, on the other hand, gradients d^u^/T, 
\d'^e\/{Te), \d'^n\/{Tn) at early times r 1 /m are large ~ 0{1), and therefore cannot be ignored. Nevertheless, 
hydrodynamics may still be applicable, provided the viscosities are unusually small r]s/seq ^ 0.1, (/sgq ^ 0.1, where 



Seq is the entropy density in local equilibrium. In this case, dissipative effects are still moderate, for example, pressure 
corrections from Navier-Stokes theory (3) 



(28) 



P \Seq S 



Heat flow effects can also be estimated based on (4) 

SKs ^ n V^(M/r) 



For RHIC energies and above, at midrapidity, the correction is rather small even for large Kg because the baryon 
density and therefore pb/T is very low. For example, in a recent ideal fluid calculation at RHIC energy [36], these 
ratios were Ub/s « 2.2 • 10~^ and Hb/T « 0.2 in order to reproduce the observed net baryon spectra. These choices 
are also supported by thermal model analyses of particle ratios which lead to hb/T ~ 0.17 [37]. 



B. Viscous equations of motion for longitudinally boost-invariant 0+lD dynamics 



At this point we specialize the equations of motion to a viscous, longitudinally boost-invariant [48] system with 
transverse translation invariance and vanishing bulk viscosity: 



e + 



risT ( 



n 

n H — 

T 

e+p 



^ n{T) = 



Ton{To) 



2r 2 \r]sT 



T 

' 3t 
" 2 



(29) 
(30) 

(31) 

(32) 



This special case is well known in the literature [9, 34, 38] as a useful approximation to the early longitudinal expansion 
stage of a heavy ion collision for observables near midrapidity « 0. Here t = Vt^ — z'^ is the Bjorken proper time, 
and the 'dot' denotes d/dr. ttl and ttt are the viscous corrections to the longitudinal and transverse pressure, i.e. 
the TTzz and = nyy components of the shear stress tensor evaluated at local rest frame[49], respectively. All 
the other components of the stress tensor arc zero due to symmetry. There is no equation for heat flow because 
the symmetries of the system — longitudinal boost-invariance, axial symmetry in the transverse plane and rj —rj 
reflection symmetry — force the heat flow to be zero everywhere. We have chosen to ignore bulk viscosity since shear 
viscosity is expected to dominate at RHIC. In the following we also concentrate on a system of masslcss particles, 
where bulk viscosity is zero in general. It is worth noticing that these equations arc identical in both Eckart and 
Landau frames, but in less restricted systems where heat flow is nonzero, Eckart and Landau frames differ. 

To simplify the discussion and to facilitate comparison with transport results, from here on we concentrate on a 
system of massless particles with only 2 — > 2 interactions. Particle number is then conserved and the equation of state 
is 



3p, 



(33) 



Now the density equation decouples entirely and we end up with two coupled equations for the equilibrium pressure 
and the viscous correction ttl. The shear stress relaxation time (12) and the shear viscosity (13) can be recast with 
the transport mean free path Xtr = 1/ {ncrtr) as 



rjs = CnTX, 



tr 5 



0.1 



(34) 



and (30)-(31) can then be written as 
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where 

The ratio of expansion and scattering timescales K controls how well ideal and/or dissipative hydrodynamics applies. 
This is essentially the inverse of the ratio of shear stress relaxation time to hydrodynamic timescales Ttt/t = 3/(2k). 
K is as well a generalization of the Knudsen number i/A, since the shortest spatial length scale is given by gradients 
in longitudinal direction L ~ XjipzUz) ~ r. It is also a measure of the shear viscosity to entropy density ratio because 
for a system in chemical equilibrium Seg = 4n and thus 



(see Sec. VE for the general case). 
Similar treatment to relativistic Navier-Stokes theory leads to 



(38) 



TTi, = -— = -— (39) 



3t 3k 
and the equation of motion 

4p 4 p 
3r 9k(t) r 

As discussed in the previous Section, the viscosities in NS and IS theories differ and, therefore, k in (40) is not identical 
to the one in (36). We will ignore the difference because in our case it is only « 5%. 
Finally, we note that in the "naive" Israel-Stewart approximation (36) changes to 

'^^ + ^^--^- (41) 



IV. BASIC OBSERVABLES 



Here we introduce the basic observables investigated in this study, and discuss their evolution based on the analytic 
Israel-Stewart and Navier-Stokes solutions of Appendix C. It is important to emphasize that our observations will 
hold only during the longitudinal expansion stage of heavy ion collisions. After some time t ~ R/cs, transverse 
expansion sets in and hydrodynamics, whether Israel-Stewart or Navier-Stokes, eventually breaks down because for 
expansion in three dimensions \tr ^ t^/ct, i.e., k ^ ct/t^ — > in the hadronic world where cross sections are bounded. 
It is interesting to note that rjs/seq ~ const would not decouple even for a three-dimensional expansion (because in 
that case T cx 1/t and thus Xtr oc rj/p cx t, while Texp = ^/{dfiU^) r, i.e., k ~ const). 

Throughout this section and the rest of the paper, the subscript '0' refers to the value of quantities at the initial 
time To (e.g., = A{to)). The most important parameters in the problem are the initial Knudsen number Kq, or 
the corresponding kq, and the initial shear stress to pressure ratio = T^Lfi/Po- 



A. Pressure anisotropy 

The magnitude of dissipative corrections can be quantified through the ratio of viscous longitudinal shear and 
equilibrium pressure 

(42) 

P 

A suitable equivalent measure is the pressure anisotropy coefficient 

which is the ratio of the transverse and longitudinal pressures px = p — tt^/S, Pl = p + ttl. In the ideal hydro limit 
the anisotropy is unity Rp ^ 1. 
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The time-evolution of the anisotropy coefficient is given by the equations of motion (35) and (36): 

" 3r (2 - ■ ^^^> 

Thus, in IS theory the pressure anisotropy is a constant of motion when the viscous stress is equal to its Navier-Stokes 
value (39), or at asymptotically late times r ^ 00. In contrast, from NS 

P 3At + 2' 

which is only constant for k(t) = const (constant cross section), or in the ideal hydro limit n ^ 00 (in which case 
Rp ^ 1). From the above it also follows that in the special case of our boost invariant scenario, if the cross section is 
constant and the shear stress starts from its Navier-Stokes value, Navier-Stokes and Israel-Stewart theory coincide. 



B. Longitudinal work 

Dissipation also affects the evolution of the equilibrium (or average) pressure. From (35), for ideal hydro evolution 
the pressure drops as p{t) oc t~^/^ due to longitudinal work. In the viscous case, the work done by the system is 
smaller because the viscous correction to the longitudinal pressure is usually negative ttl < 0. Therefore, pressure 
decreases slower than in ideal hydro, and deviations from the ideal evolution, such as the ratio 

p(t) T(t) 

~ (for conserved particle number) (46) 



Pideal (t ) Tide 

can be used to quantify dissipativc effects. 

Studies in the past [27, 28] have analyzed a closely related quantity, the transverse energy per unit rapidity, dET/drj. 
This is simply a combination of the pressure anisotropy and deviation from ideal pressure 

dEr _ 3nT dN / 5 A _ 37rro dAT /tq \ - V3 p(t) 3[7 + Rpir)] 
dr] 4 dr] \ 16V 4 dr]\T) Pideai {t) 8[2 + Rp{t)] ^ ' 

provided the quadratic ansatz (11) is applicable (see Appendix D 1). 

We can make a few generic observations based on the analytic Isracl-Stcwart and Navier-Stokes results (C3)-(C4), 
(C8), (C21)-(C22), and (C29) from Appendix C. For a constant cross section, p/pideai grows without bound - 
dissipative corrections keep accumulating forever. The influence of the initial shear stress, or equivalently shear stress 
to pressure ratio ^0 = C('''o), is of C'(^o/'«o) and thus vanishes in the large kq limit. At late times r ^ tq, for Kq > 2 
and not too large initial shear stress to pressure ratio |fo| <C 2kq 



Pideal / „=const V "^0 / 9ko V / 3^0 3Kq 2ko 

i.e., for T « lOro and Kq — 2 the accumulated pressure increase is p/pideal ~ 1-3, while p/pideal ~ 1-15 if Kq = 5. For 
a scale invariant system with rjs/seq ~ const, on the other hand, dissipative effects are more moderate for the same 
Kq and at late times approach a finite upper bound 



Pideal / rj/s^const 



A_ M ^ i + (49) 

3ko 2/to / 3ko 2k;o 



This is because scale invariant systems turn more and more ideal hydrodynamic as they evolve (as long as their 
expansion is only longitudinal). For the same Ko = 2 and 5 with ~ 0, the bounds are modest, p/pideal ^ 1-25 and 
< 1.1, respectively. 



C. Entropy 



Another quantitative measure of the importance of dissipative effects is entropy production. Here we consider an 
ultra-relativistic system (thus 11 = and (32 = 3/(4p)) with 2^2 interactions, ID Bjorken boost invariance, and 
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transverse translational, axial, and rj — > —r] reflectional symmetries (imply =0). Therefore, the entropy current 
(14) simplifies to 

S^ = su^, s = s,,-^ (50) 

where 

71/ Li 

.e, = n(4-x), x-ln^ = ^ (51) 



and 



ne,(r) = (52) 



is the particle density in chemical equilibrium at temperature T for massless particles of degeneracy g in the Boltzmann 
limit. Dissipative contributions in the entropy density s are negative, in accordance with the fundamental principle 
of maximal entropy in equilibrium. 

The equations of motion (35)-(36) imply an entropy production rate 



d,S^ = ^drirs) = > . (53) 

r 4t 



Equivalently, the entropy per unit rapidity 



never decreases 



^^tAts (54) 
dr] 



Here At is the transverse area of the system, and in the last step we substituted the rapidity density dN/drj = tAtu. 
Equation (54) is a special case of a generalized conservation law (B7) applied to the entropy current S^^ 
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dx^ d^Sf" = dr(j J dxl S^^^ - J dxl S^^ . (56) 



Analogous relations can be obtained for the energy, momentum, and charge density. In 0+lD these are quite trivial 
- they respectively reproduce (35), give identically zero, and dN/dri = const. In higher dimensions, however, the 
generalized conservation laws present important constraints that any solution must satisfy at all times and, therefore, 
they are ideal for testing the accuracy of numerical solutions at each time step (see Appendix B). 

Only the complete set of Israel- Stewart equations gives the correct rate of entropy production. The 'naive' approxi- 
mation does not guarantee a monotonically increasing entropy 

\naive IS _ SkH. 2 f ^ + 4 



(a^S'^)™ = (^1 - , (57) 

unless K is sufficiently large and, away from equilibrium, it underpredicts for a given ^ the entropy production rate [50] 
(since ^ < — 1 is unphysical). In contrast, the second law of thermodynamics does hold for Navier-Stokes 

id.ST' = ^e^s > . (58) 

The NS result is the same as (53) but with the shear stress restricted to its Navier-Stokes value. We note that in 
Israel- Stewart theory the naive entropy per unit rapidity, defined using the equilibrium entropy density 



dS' 

dr} 



= Seg tAt (59) 



does not increase monotonically. Rather, it increases (decreases) for negative (positive) ttl- 
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Based on the analytic Israel-Stewart and Navier-Stokes results in Appendix C, we can outline general expectations 
for the entropy evolution. For a constant cross section by late times r ^ Tq the entropy increase relative to the initial 
entropy is logarithmic with time 



(dS/drj) 



{dS/d'n)o 



1 



(T=const 



4-Xo 



3 In — 

Pideal 16 



1 



4-Xo 



T 3 16 3^0 



To 



(60) 



where we considered initial conditions not too far from local equilibrium. E.g., by r « IOtq with Kq = 2 and chemical 

equilibrium initial conditions w 20% entropy is produced, while w 10% with Kq = 5. For a scale invariant system 
with rjs/seq = const, on the other hand, entropy production is slower for the same Kq and saturates at late times 



[dS/dr]) 



{dS/dr])o 



1 



1 



7]/ s^const 



4 - xo Ko 



1 



(?) 



2/3 



3Co 

4 



1 



XO Ko 



1 



3^ 
4 



Vs 



TqTo s 



eg 



1 



3^0 
4 



(61) 



For the same Kq = 2 and 5 (and ~ 0), the entropy increase by r = IOtq is smaller, w 15 and « 6%, respectively. 
Based on this simple analytic formula for entropy production, we also confirm the results of Ref. [38] , which considered 
IS hydrodynamics with a unique initial condition ~ — 16/(9Toro) x rjs/seq where To w 0.39 GeV x (0.14 fm/ro)^/^ 
and To was varied between 0.5 and 1.5 fm. 



V. REGION OF VALIDITY FOR DISSIPATIVE HYDRODYNAMICS 



Here we determine the region of validity of dissipative hydrodynamics via comparison to full nonequilibrium two- 
body transport theory [26-29]. We consider two scenarios: Scenario I with a constant cross section, which is least 
favorable for hydrodynamics; and Scenario II with a growing a oc r^^^, which is the most optimistic for applicability 
of hydrodynamics and is very close to rjs/seq = const as we show in Appendix C. In the same Appendix we also study 
a scenario with a cx l/T^ that turns out to be close to Scenario II but with stronger dissipative effects, and discuss 
analytic Navier-Stokes and (approximate) Israel-Stewart solutions. 

Due to scalings of the equations of motion, the results presented here are rather general. Equations (35)-(36) are 
invariant under rescaling of time, and/or joint rescaling of the pressures p and ttl, provided the dimensionless k 
depends only on p, tt^, t/tq and no additional scales (all solutions studied here satisfy this condition). The same 
scalings are exhibited by the transport [28]. For a physically reasonable po > 0, it is therefore convenient to consider 
dimensionless pressure variables p{t) = p{t)/po and i^l{t)/po, for which the solutions only depend on f = t/tq, 
Kq = Kq/C and the initial condition = t^l.q/po- 

Unless stated otherwise, we initialize the transport based on the quadratic form (11). In our case of an ultrarela- 
tivistic system (e = 3p) in the Boltzmann limit with vanishing bulk pressure and heat flow 

D^ = 0^ C^^ = _ ^ ^^(^ = 0,p) = ^^^;/i, (62) 



where p± = ^P^+Py is the transverse momentum. We ensure nonnegativity of the phase space distribution via a 
9-function 

/(,? = 0, p, r = ro) = 1^ ^ 1^ [1 + </-G(r?, p)] 6(1 + <t>G{r,, p)) (63) 

where At is the transverse area of the system (with the elimination of negative phase space contributions, a normal- 
ization factor F(^) < 1 is needed to set a given dN /dif). The cutoff does not affect the general scalings of transport 
solutions but does influence the initial pressure anisotropy (for example, values Rp = 0.3 and 1.75 set based on (62) 
change to Rp « 0.476 and 1.693 when the cutoff is applied). Therefore, we initialize hydrodynamics with a shear 
stress TTi that gives the same initial pressure anisotropy as the transport. 

The transport solutions were obtained using the MFC algorithm [39], which employs the particle subdivision tech- 
nique to maintain covariance[26, 28]. Transverse translational invariance was maintained in the calculation through 
periodic boundary conditions in the two transverse directions. A longitudinal boost invariant system was initialized 
in a coordinate rapidity interval — 5 < 77 < 5, and observables were calculated via averaging over — 2 < 77 < 2 with 
proper Lorentz boosts of local quantities to 77 = 0. 
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FIG. 1: Time evolution of pressure anisotropy Rp = pl/pt from covariant transport (solid lines with symbols) and Israel-Stewart 
dissipative hydrodynamics (solid lines) as a function of K = r/Atr(r), from local equilibrium initial conditions nL{To) = 0. 
Results for Navier-Stokes (dotted lines) and free streaming (dotted line with circles) are also shown. Left: a = const scenario, 
for which the curves are labeled by K{r) = const = Ko ~ 1, 2, 3, 6.67, and 20. For K = 1, the Navier-Stokes result is 
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Ko = K{to) = 1, 2, 3, and 6.67. 



Right: 



2/3 



scenario, for which rjs/seq ~ const and the curves are labeled by the initial 



A. Pressure anisotropy 

Figure 1 shows the pressure anisotropy pl/pt evolution as a function of the rescaled proper time f = t/tq from 
the transport (solid lines with symbols) and Israel-Stewart hydrodynamics (solid lines without symbols) with local 
equilibrium initial condition. The left panel shows calculations for the a = const scenario. For Kq = 1, the anisotropy 
from IS hydro starts to fall rapidly below the transport above t > 2to and it is a factor ~ 5 smaller by late r ^ lOro. 
Clearly, the system cannot stay near equilibrium when the rate of scatterings equals the expansion rate. With 
increasing Kq, the undershoot becomes smaller and gradually vanishes as Kq — > oo. The difference is only ~ 10% 
already at Kq = 3, and is rather small by Kq « 7. 

The right panel shows the same but for the growing cross section scenario with rjg/ S(.q ~ const. The situation of 
course improves because in this case K increases with time. For Kq = 1, IS hydro undershoots the pressure anisotropy 
from the transport only by ~ 20% and the differences vanish at late times (since in this case both theories converge 
to i?p = 1 as r —> oo). ~ 10% accuracy is achieved already for Kq — 2, while for Kq = 3, IS hydro is accurate to a 
few percent. 

Moreover, the above findings hold for a wide range of the initial conditions, including large initial pressure 
anisotropies, as shown in Figures 2 and 3. These figures are for the same calculation but with i?p(ro) = 0.476 
and 1.693, respectively (which correspond to = —0.423 and 0.375). We emphasize that the results hold only if 
noncquilibrium corrections are close to the form (11) suggested by Grad. For such class of initial conditions, however, 
we find that Israel-Stewart hydrodynamics can well approximate the transport 10% accuracy) provided Kq > 3, 
even for the most pessimistic constant cross section scenario. If rjs/seq — const, only Kq > 2 is needed. We stress 
that in either case, there is no need for the initial conditions to be near the Navier-Stokes limit. 

This is quite remarkable because from Figs. 1-3 it is clear that already the early evolution differs between IS 
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hydrodynamics and transport. E.g., for an equilibrium initial condition (^(tq) = 0), IS hydrodynamics (44) gives 

i?f(r) = l-^^^+0((r-ror) (64) 
for any initial value and evolution scenario for k. From covariant transport, on the other hand (see Appendix D 2) 

i?*™-P(r) = 1 - fc_Z£l + 0{{t - Tof) . (65) 

5ro 

I.e., pressure anisotropy develops, universally, 20% faster from the transport than from IS hydrodynamics (if the 
evolution starts from equilibrium). 

This illustrates a limitation of the hydrodynamic description of transport solutions. Similar discrepancies were 
observed in [8] in the early evolution of differential elliptic flow V2(jpt)- Remarkably, in our case, though the trans- 
port develops deviations from equilibrium faster, its rate of departure slows down quicker, which at intermediate 
times results in smaller accumulated dissipative corrections to the pressure anisotropy than from IS hydrodynamics. 
Eventually, the hydrodynamic evolution "catches up" to the transport, except for low K <3m the a = const scenario. 

Figures 1-3 also show the Navier-Stokes approximation (dotted lines without symbols) for each of the Israel-Stewart 
results. By late times, the Navier-Stokes and Israel-Stewart solutions converge for both cross section scenarios, 
independently of the initial pressure anisotropy (for a = const and Kq = 1, the NS anisotropy is negative and 
therefore not visible in the plots). However, the applicability of Navier-Stokes theory at early times depends, besides 
the value of Kq, strongly on how far the initial shear stress is from its Navier-Stokes value (39). Navier-Stokes assumes 
that shear stress, and therefore the pressure anisotropy, relaxes immediately, but relaxation happens over a finite time. 
The approach toward the Navier-Stokes limit is governed by = 3t/ (2k), therefore Navier-Stokes becomes applicable 
only after some time At ~ |i?o — Rns\ tq/k. Note that the initial slope of the R{t) curves does not always reflect r^r 
directly because it is given by the initial derivative of ^ 

R{r) - |^(t) = --^(C - iNs) + 0{1)^ (66) 

where wc combined (26), (35), the observations that Y ^ 0{\)/t and Z = 0, and assumed is small. For local 
equilibrium initial conditions the slope of R{t) is therefore ~ C(l)CArs/T7r ~ 0{1)/t, independently of Kq (cf. 
Figure 1 and also (64)). For initial shear stresses far away from the Navier-Stokes limit, on the other hand, the slope 
~ 0{1)^/tt^ cx k steepens with increasing K as seen in Figures 2 and 3. 

The inaccurate description of early shear stress evolution in Navier-Stokes has a cumulative effect on the evolution 
of thermodynamic quantities, such as the pressure and the entropy, as we show in the next two Sections. Of course, 
the errors are proportional to ratio of the time the system spends away from the NS limit and the hydrodynamic 
timescale, i.e., At/tq ~ 1/k. 

B. Pressure evolution 

Now we turn to the evolution of the (average) pressure. In ideal hydrodynamics {Kq oo) the pressure drops rapidly 
with time pid a t~^/^. Therefore it is more convenient to study dissipative effects relative to ideal hydrodynamics 
through the ratio p{t) / pidir) . 

Figure 4 shows the pressure relative to that in ideal hydrodynamics as a function of the rescaled proper time f = t/tq 
from the transport (solid lines with symbols) and Israel- Stewart hydrodynamics (solid lines without symbols) with 
local equilibriTim initial condition. The left panel shows calculations for the <j = const scenario. For all Kq values, 
the evolution starts out the same between IS hydro and transport but then the hydro starts to accumulate deviations 
because it follows the shear stress evolution only approximately. For Kq = 1, IS hydro maintains 10% accuracy in 
the magnitude of dissipative corrections (i.e., p/pid — 1) only up to r « Atq. As Kq increases, the situation improves 
gradually, for = 3, 10% accuracy holds up to r « IOtq, and by /To w 7 the hydro stays within a few percent of 
the transport even until r — 2Qtq. 

The right panel shows the same but for the growing cross section scenario with rjs/seq « const. This scenario is 
more favorable for the hydrodynamic approximation because K ~ r^^^ grows with time. For Kq = 1, the error in the 
dissipative correction [p/pid — 1) is less than 10% up to t w Stq, and already for Kq = 2 IS hydro is accurate to within 
better than 10% throughout the whole range t < 20to studied. The pressure evolution results therefore reinforce the 
regions of validity found in the previous Section {Kq > 3 for = const, and Kq>2 for r]s/Sf,q ~ const) 

Clearly, the region of applicability for Navier-Stokes is more limited (Figure 4, dotted lines without symbols). For 
low Kq, it overestimates the pressure corrections not only at late times but also at early t ~ few x tq. Kq « 7 is 
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FIG. 4: Same as Fig. 1 except for the time evolution of the (average) pressure. The pressure is plotted normalized to the 
pressure Pideai{T) — Po{to/tY^''^ in ideal hydrodynamics. 

barely sufficient for 10% accuracy in viscous corrections for rjs/ Seq ~ const, but it is not enough in case of cr = const. 
Based on the trends with increasing Kq, we estimate that Kq > 9 — 10 is needed for Navier-Stokes with a = const to 
deviate less than 10% from the viscous effects calculated with the transport. Therefore, for local equilibrium initial 
conditions, Navier-Stokes theory becomes applicable at about three times shorter mean free paths, or equivalently three 
times larger longitudinal proper time t (i.e., three times slower longitudinal expansion) , than Israel- Stewart theory. 



C. Entropy 



Now we proceed with results on entropy production. In transport theory, the entropy current is defined as 

d?p 



\n{^^f{x,p) 
9 



(67) 



where g is the number of internal degrees of freedom. This nonlinear function of the phasespace density / is cumber- 
some to evaluate with the MPC code, and therefore we here opt for an approximate result based on the truncated 
Israel- Stewart expression (50), evaluated using the pressure and shear stress from the transport. This includes dissi- 
pative corrections to the entropy up to quadratic order in (/). 

In the most dissipative a — const scenario with Kq ~ 1, there is about 30% additional entropy produced by late 
times t/tq ^10—20 as can be seen in Figure 5 (left plot). For rjs/seq ~ const (right plot), the same Kq — 1 yields 
only about 20% extra entropy. With increasing Kq entropy generation gradually weakens and by Kq ~ 7 it is only 10 
and 5 %, respectively. 

The Israel- Stewart results are within 15% of the approximate transport results already for Kq = 1, and about 10% 
accuracy in the calculated dissipative effect is achieved for Kq > 3 (for a = const) and Kq > 2 (for rjs/seq ~ const). 
In contrast, the Navier-Stokes strongly overpredicts the entropy, unless Kq exceeds about 6 for cr = const or « 3 for 
Vs/seq ~ const. The bounds for 10% accuracy are in agreement with those found previously in Sec. VB. 



D. Limitations of the 'naive' Israel-Stewart approximation 



Now we discuss the applicability of the 'naive' Israel-Stewart equations. Figure 6 compares the pressure evolution 
in complete Israel- Stewart theory to that in the naive approximation, for local equilibrium initial conditions (^o = 0), 
as a function of the rescaled proper time t/tq. Clearly, the naive result overshoots the pressure both for the constant 
cross section scenario and for rjs/seq ~ const, unless Kq is large. This confirms expectations based on the analytic 
solutions in App. C. Though the 'naive' theory converges to the correct result at large enough Kq ^ 7—20, comparison 
with Fig. 4 tells that it is even less accurate than Navier-Stokes theory. 
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The reason for the large errors is that away from local equihbrium the 'naive' approach drives the shear stress more 
negative (compare (36) and (41), and note that typically tt]^ < 0). This is demonstrated in Fig. 7 where we plot the 
pressure anisotropy Rp, which is a monotonic function of ^ = ttl/p- For cr = const, we find that the naive approach 
saturates the anisotropy at a lower value than the complete theory, confirming analytic expectations in App. CI. 
For rjs/seq ~ const, the system does approach ideal hydrodynamic behavior eventually, however that occurs on a 
much longer timescale than from complete Israel-Stewart theory. This is in agreement with expectation based on the 
analytic solutions (C34)-(C36). 

The pressure anisotropy results further reinforce our conclusion that the 'naive' Israel-Stewart approximation is 
poorer than Navier-Stokes (cf. Fig. 1). In heavy-ion collisions, gradients are large, at least initially, and therefore 
cannot be ignored even if dissipativc corrections (e.g., t^l/p) are small. 
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E. Implications for heavy-ion physics 

Having determined the region of validity (defined as 10% accuracy in dissipative effects) for Israel-Stewart and 
Navier-Stokes hydrodynamics in terms of the initial ratio of the expansion and scattering timescales Kq = To/Xtr,Q 



(cr ~ const) 

iVs/Seq ~ const) 



we now turn to implications for heavy-ion collisions. From (37), (51), and (52), 
TqTq sq ^ 1 f To \ f To \ fl/i^n) 



Xo Vsfl 



15.9 X 



l-Xo/4 Vl GeV 



\1 fm/ 



Ko w 0.8ko 



(68) 
(69) 



(70) 



Therefore, we can place an upper limit on the (initial) shear viscosity for which IS or NS reproduces with better than 
10% accuracy the viscous corrections to basic observables such as pressure and entropy: 



Seq.O 



< O.SToTo 



IS 



Seq,0 

47r?7s 



< 0.25roTo 



NS 



IS 



< 0.4Toro 



(cr ~ const) 

{Vs/Seq ~ const) 



(71) 
(72) 



NS 



where we assumed chemical equilibrium initial conditions {xo = 0). If the shear viscosity of dense quark-gluon matter 
is bounded from below by A-Krjs/seq > 1, as has been conjectured recently, then the situation for Israel-Stewart is 
close to marginal. For rjs/seq = l/(47r), typical parton transport initial conditions (Tq = 0.7 GeV, tq — 0.1 fm) 
translate into Kq < 1, for which IS is not applicable for either of Scenario I or II, while for typical hydrodynamic 
initial conditions {To ^ 0.38 GeV, tq = 0.6 fm) we have Kq < 3, sufficient for both scenarios (barely for a = const). 

On the other hand, Navier-Stokes may be marginally applicable only if rjs/seq ^ 0.5/(47r) throughout the whole 
evolution, at least based on this 0+lD study, where acausal artifacts and instabilities do not arise. We emphasize 
that the bound quoted here is for initial conditions close to local equilibrium. The accuracy of the Navier-Stokes 
approximation strongly depends on how far the initial shear stress is from the Navier-Stokes value. If the evolution 
starts out near the Navier-Stokes limit, we expect Navier-Stokes to be accurate up to higher viscosities. 

Within the region of applicability of Israel-Stewart, dissipative corrections to the average pressure and the entropy 
are modest and stay below ^ 20% even up to late times r < IOtq. This may serve as a useful "rule of thumb" 
applicability condition for hydrodynamics: if dissipative corrections to average pressure and the entropy calculated 
from hydrodynamics are significantly larger than 20%, the validity of hydrodynamics is questionable. 

The above findings reinforce a recent calculation[18] in 2-l-lD that found good agreement between IS hydrodynamics 
and 2 — + 2 transport, for conditions expected in Au + Au at y/SNN ^ 200 GeV/nucleon at RHIC, in case of a small 
shear viscosity to entropy density ratio rjs/seq ~ l/(47r) (on average). The same study also found good agreement 
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between the two theories for a large constant transport cross section atr ~ 13 mb. That is also in line with our results 
here because it corresponds to 47r?7s/sgq(To) ~ 0.25 in the center of the collision zone, i.e., initially rjg/ s^q < l/(47r) in 
most of the system. 

Finally we note that the applicability of the hydrodynamic approach on very short time and length scales is 
another important question. In typical real-life problems TqTo ^ 1 because the hydrodynamic expansion timescale 
r is by orders of magnitude larger than the quantum (energy) timescale 1/T. This also leaves ample room to make 
hydrodynamics applicable (kq ^ 1) even for appreciable viscosities. In the heavy-ion case, however, the two timescales 
are comparable Tqto ^ and therefore a macroscopic treatment may be marginal. 

VI. CONCLUSIONS 

Based on comparison to covariant transport theory, we explore the region of validity of Israel-Stewart and Navier- 
Stokes hydrodynamics in heavy-ion physics applications. We follow the evolution of the average pressure, pressure 
anisotropy, and entropy for a massless ideal gas in 0+lD longitudinally expanding Bjorken geometry. Binary 2 — > 
2 interactions are considered for two main scenarios, a fixed cross section a = const (Scenario I, pessimistic for 
hydrodynamics) and a scale invariant system with rja/sgq « const (Scenario II, optimistic for hydrodynamics). 

We find (Sec. V) that dissipative effects calculated from Israel-Stewart hydrodynamics reproduce those from the 
transport to within 10%, provided initially the expansion timescale is three (for Scenario I) or two (for Scenario II) 
times larger than the transport mean free path, i.e., the initial Knudsen number Kq = To/Xtr,Q 3 or 2. When this 
criterion is fulfilled, Israel-Stewart is accurate even if initial pressure anisotropies are large Pl/pt ~ 0.4 — 1.7 - there 
is no need to start near the Navier-Stokes limit. On the other hand, same accuracy from Navier-Stokes requires three 
times larger Kq, if the expansion starts from local thermal equilibrium (unlike for Israel-Stewart, the applicability 
of Navier-Stokes depends strongly on how far the initial shear stress is from its Navier-Stokes value). We emphasize 
that these findings apply only when initial viscous corrections are of the quadratic form suggested by Grad (11). 

These results imply that (Sec. VE), for typical heavy ion initial conditions at RHIC energies, Israel-Stewart hydro- 
dynamics is accurate up to rjg/seq < 1.5/(47r), while for Navier-Stokes rjs/seq < 0.5/(47r) is needed. This is supported 
by a recent 2+lD calculation[18] that finds good agreement between Israel-Stewart and transport for r/s/seg ~ I/C^ti"), 
and also for a large atr ~ 13 mb. 

In addition, we test the accuracy of the naive Israel-Stewart approximation (Sec. VD) that neglects products of 
gradients and dissipative quantities in the equations of motion, and find that it has an even more limited applicability 
than Navier-Stokes. 

We also compare in detail (App. C) Israel-Stewart and Navier-Stokes solutions in 0-l-lD for four scenarios, a = const, 
(T oc 1/T^, (J oc r^/''* and rjg/seq = const, and find that results for the latter two arc almost identical, even at low initial 
Knudsen numbers Kq ~ 1. Moreover, we obtain analytic Israel-Stewart and Navier-Stokes solutions in 0-l-lD, which 
are useful for quick estimates (Sees. IV B and IV C) and to test numerical solution techniques. We also derive additional 
tests (App. B) based on generalized conservation laws for conserved currents, energy-momentum, and entropy, which 
can be utilized to verify the accuracy of numerical Israel- Stewart solvers in 1-1-1, 2-1-1, and 3-1-1 dimensions. 

Finally we emphasize that the current study is limited to a massless ideal gas with particle number conserving 
interactions in 0-l-lD Bjorken geometry. The infiuence of the transverse expansion will be quantified in a future paper 
(requires at minimum a 1-l-lD approach). It will be also important to check how the results depend on the equation 
of state and the presence of particle non-conserving processes, such as radiative 2^3. For a nonconformal equation 
of state, bulk viscosity may become important[41, 42]. Ideally, one should also test the accuracy of the hydrodynamic 
approximation for nonequilibrium theories other than covariant transport. 
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APPENDIX A: ORIGIN OF oo, oi, o^, oi IN THE ISRAEL-STEWART EQUATIONS OF MOTION 

The equations of motion (23)-(25) reproduce the entropy production rate (16) only approximately, up to typically 
small quartic and higher- order corrections in dissipative quantities. With ai = = a'-, a contribution 

Uq^'TV^iao/T) - q^n^^^TV ^{a^/T) (Al) 

would be missing from Td^S^^ in (16). These terms are bilinear in the dissipative quantities and, therefore, can be 
split arbitrarily between the bulk pressure and heat, and heat and shear equations of motion. I.e., with 

TW{ai/T) = af + a'/" (A2) 

(16) is identically satisfied but contributions to the equations of motion depend on aj 

PoDH = (...)+ a;/^g^ (A3) 
piDq^ = (...)+ A>o"n-« (A4) 

/^a^'TT"" = (...)- a'/"g''> (A5) 

Only components orthogonal to contribute but apart from that constraint Oq and are arbitrary functions 

of the hydrodynamic fields and their derivatives, and potentially new scalar functions {a\''\e,n)} characterizing an 
isotropic matter. However, ignoring the dependence on the dissipative quantities is consistent with the truncation of 
the entropy current (14) at quadratic order in those. Moreover, for small deviations from equihbrium one may seek 
to include only the leading contributions coming from first derivatives of the ideal hydrodynamic fields, i.e., 

< = al'^Du-^ + af )rV^^ + V^^ (A6) 

where wc chose 1/T and n/T as the two independent variables instead of e and n. But the three terms arc not 
independent - energy- momentum conservation (1) and the Gibbs-Duham relation sdT = dP — ndjj provide one 
constraint 

^A'^'DUa + V";^ = ^^V"^ , (A7) 

T T e + p T ' ^ ' 

and W^{ii/T) may be ignored, at least parametrically, because it is proportional to the heat flow (4) in the first-order 
(Navier-Stokes) theory. Therefore, to leading accuracy only one scalar function enters and we can write 



Analogous arguments give 



= -ai{e,n)Du^ . (A8) 



V ia,/i ; ~ ^^^^^^ V ^ ~ g^^^^^ V nu„ (Ayj 



from which (20) follows. 

We plan to revisit the above approximations in a future study. In any case, they do not influence our 0-l-lD 
calculations here because the terms do not play a role (heat flow vanishes by symmetry). 



APPENDIX B: GENERALIZED CONSERVATION LAWS 



Here we present general relations of the form 

^=B(r) (Bl) 

that can be used to test the accuracy of numerical dissipative hydrodynamics solutions in any dimensions. A and B 
only depend on the hydrodynamic fields at the given r. Evaluating them at each time step, one can either numerically 
differentiate A{t) or integrate B{t) and check how accurately the solutions satisfy (Bl). 
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Consider a four- divergence d^A^^{x) (in regular Minkowski coordinates). Integration over a four- volume V4 gives 

jd\d^A^'{x)= j da^,{x)A^^{x) (B2) 

Vi (T{Vi) 

where (7(14) is the 3D boundary ("surface") of V4. Now take the special case of a Bjorken " box" V4 = At x Ary x At 
with an infinite transverse area At —f 00 but infinitesimal proper time and finite coordinate rapidity extensions 
At — * 0, At; = t]2 — Assuming A'^{x) drops faster than l/x|, at large |xt|, we can neglect surface terms at 
\xt\ — > 00 and keep only contributions on r = const and r] = const hypersurfaces: 



/ 



drrdr? da;|a^A^(x') 



da(-) - J da(-) + J da^") - J d^'') 

<7(r+dr) o-(r) cr{ri2) "'('Ji) 



Af'ix) (B3) 



where the surface normals are 



dcr^^^ = rdxTdrjUg , '^'^(ri) ^ — drdxy Ug , with = (chry, 0, shr;) , = (sh?7, 0, chr;) (B4) 

and we used d'^a; = drTdrjdx^. Here, is the longitudinal Bjorken flow velocity, while U3 is its orthonormal 
counterpart in the t — z plane. Note that the actual flow velocity does not need to be us. Dividing by At and taking 
the limit we arrive at 

T J dr?d4 a^A^ ^dr(^T J d77d4 u%A^^ - J d4 {AM - Mm)) , (B5) 
which is a generalized conservation law for the quantity 

■^=T j drydxK^A* • (B6) 

If dfj.A'^ = 0, and the surface term 1*3(^^(771) — ^^(772)) vanishes, we have A{t) = const. 

In a boost-invariant calculation the longitudinal extension of the system is formally infinite and thus a generalized 
conservation law for a quantity per unit rapidity is more practical. It can be obtained in a similar fashion if one 
divides by A77 and takes the limit At; 0. The result is 

T J d4 ^^A^' = dr^-dr,J da;| u^A^ , (B7) 



^ = T / d4u^^^ . (B8) 



where 

dr] 

Again, if d^A^^ = and the ?7-derivative term vanishes, we have dA/drj = const. 

1. Charge / Particle number 

We first apply Eq. (B5) to a conserved current in Eckart frame: A^'^ = ngqU^, where m'^ = 'y{ch9,v eii,sh9) is the 
flow four- velocity and 9 is the flow rapidity. Now u'^U/j, = 7ch(ry — 9) and u'^Un = 7sh(ry — 0). If the rapidity interval 
is so large that A'' (r/i) = N^{ri2) = 0, or the system is boost invariant, r] = 9, the surface terms are zero and we get 
a simple conservation law 

N = T J dr] d4 7 cli(?? — &) = const . (B9) 
In a boost-invariant case, the coordinate rapidity integral is trivial and we get 

d47'^ = const . (BIO) 



— -T /( 

dij J 
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2. Entropy 

Second, we apply Eq. (B5) to the entropy current (14) and its divergence (16). If S'''(r/i) = S''{t]2) = 0, we get 

a,s = ./d.,d,^(|-i^ + ^)>o, (BID 

where the entropy of the system is 

S = T J dr/da;|u^S'^ , (B12) 

and the last inequality follows from the general properties q'^Qfj. < and tt^^'^tt^^ > 0. 

For longitudinally boost invariant dynamics, it is more natural to follow entropy per unit rapidity: 



drj J J^^' \dr]J J ' \CT k^T"^ 2risT 

3. Energy-momentum 

Finally we derive the conservation equation corresponding to energy-momentum conservation dpT^'^ = 0. Contrac- 
tion of the energy-momentum tensor with gives the conservation of energy. In a case where the entire system is 
within the interval [?7i , , 

drE = dr(^T J dr/dx^u^T^^u^^ = . (B14) 

Contraction with = (0, br, 0) gives the change in transverse radial momentum. Substituting 

d^^iT^^^UR,^) = + T'^'d^UR,^ (B15) 

into Eq. (B5) results in 

drMr = dr(T J df] dx^ u^Ti^^u^r^ = t J dr? da;| T'^^a^u^ . (B16) 

To be more specific, we also show as an example a boost-invariant, cylindrically symmetric case. In Landau frame 
T^^- = (^e+p + U)uV -{p + Il)g'"' + (-^2 - ^sX^i + fr2« + n^u^u^s > (B17) 

where Ui is the orthonormal counterpart of the flow velocity in the time-radial plane, while U2 and U3 are orthonormal 
counterparts of these in the axial and beam (rapidity) direction 

u'^ = ^{chri,v BR^shrj) , = j{v chri,eR,v shrj) , = (0, e^, 0) , Ug = (shr;, 0, chr;) . (B18) 

These vectors are normalized to = 1, = w| = m| = — 1. The viscous pressure tensor components in the fluid 
rest frame are ttJ^^ = diag{0,—Tf2 — 7^3,^2,7^3). It is important to notice that the surface terms in Eq. (B5) or the 
ry-derivative term in Eq. (B7) are now nonzero. Contraction by u'^ as above and substitution into Eq. (B7) gives the 
evolution of the energy per unit rapidity: 

dr (^^^ =dr(^T j dxl T°°(ry=0)^ = - j dxlip + n + n^) . (B19) 

Contraction by gives the evolution of transverse radial momentum per unit rapidity: 

<^Mr \ ( _ f ^^2 ^01. n /■^.2 P + n + #2 



where we have used the relations 

< V< = -T°\v=0, ^=0) , d^R = --^«2,^< . (B21) 

The above results reflect general expectations. Particle number, per imit rapidity dA'^/d?7, is strictly conserved in 
both the ideal and the dissipative case. Entropy per unit rapidity dS/drj is conserved for an ideal fluid but increases 
if there is dissipation. In both cases, the energy per unit rapidity dE/dr} decreases due to longitudinal work, while 
the radial momentum per unit rapidity dM^/d?] increases due to build-up of radial flow, as long as the system stays 
near equilibrium (i.e., the total pressure is dominated by the ideal part). 
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APPENDIX C: VISCOUS SOLUTIONS FOR VARIOUS CROSS-SECTION SCENARIOS 

In this Section we analyze viscous Israel-Stewart and Navier-Stokes solutions for four different types of cross section: 
constant, a oc a oc r^/^, and Jjs/Sgg = const. For convenience, we will often use normalized quantities 



(CI) 



We will show that for typical observables of interest (average pressure, pressure anisotropy, entropy, shear viscosity 
to entropy ratio), rjs/seq = const dynamics is well approximated by cr oc t^^^ already for Ka = 1. 

In analytic considerations, it will be often convenient to drop the 7r|^ term in the equations of motion (35)-(36), 
which is a good approximation for Ittl] "C p, i.e., the general region of validity of viscous hydrodynamics. This should 
not be confused with the "naive" Israel-Stewart approximation, which also ignores the 4/3 factor in (36). For the 
a (X T^/^ and a = const scenarios wo obtain this way accurate approximate analytic Israel-Stewart solutions. We also 
derive analytic Navier-Stokes solutions for a = const, a oc r^/^ and cr oc l/T^. 



1. Solutions for ultra-relativistic gas with constant 2 — > 2 cross section 



For a constant cross section, 



Atr(r) oc T 



K{t) = 



Kq = const . 



Atr(To) 

If we ignore 7r| term, the linear equations of motion (35)-(36) can be solved in a straightforward manner: 

TTL ,0 



4 "■0 

7rz,(f) = f-3- — 



4 '^0 

P(f) = f~3~"3" 



2 T,(f)-^U.,o + f^)T_(f) 



^T+{f) + ^{KoPo-7rL,o)T.{f) 



where 



Ko = 



C 



D 



^0 ' 



T±{x) = x'^'^±x 



-D/3 



V{tq) = PO , -""/.(to) = TTLfl 



(C2) 

(C3) 
(C4) 

(C5) 



For a practical approximate formula for the pressure evolution, sec (48). 

In the ideal hydrodynamic {rjs 0, or equivalently kq oo) limit we recover 



ttl{t > to) = , p{t) = po (J^^ 



4/3 



(C6) 



At late times the pressure anisotropy, irrespectively of its initial value Rpfl, approaches a constant determined solely 
by the parameter kq 



Rp{T oo) 



12ko - 10 
9D + 3ko + 14 



< 1 



(C7) 



For a finite kq, the final anisotropy is below unity. 

Therefore, with a constant cross section, the late-time behavior of the system does not become ideal hydrodynamic 
but instead the Navier-Stokes limit applies (cf. (44) and (45)). Indeed, for large kq, (C3)-(C4) reproduce the NS 
solution 



p^^{r)=po{^) 



4/3-4/(9ko) 



4p^^(r) 
3ko 



(C8) 



and the final IS and NS anisotropics (C7) and (45) agree, i?oo = 1 - 2/ko + 4/(3k§) + 0{1/kI). 

Because i?oo is a monotonically increasing function of kq, the final pressure anisotropy is a measure of the viscosity. 
Inverting (C7), 



5 + URo 



R 



6-3i?o 



3i?^ 



(C9) 
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i.e., near equilibrium {kq » 1) 



n(T) K 2 \''o/ ' 3 9 TT-o To To 

where in the last step we approximated the temperature evolution using the leading NS term (C8). It is natural to 
measure viscosity relative to the density, which up to a factor (4 — x) is the same as rjs/seq- 

The exact analytic solutions to the 'naive' Israel-Stewart equations are analogous to (C4)-(C3) but involve different 
powers of f 



2 Ato ^ V'^o-^^o + f 

The late time behavior is governed by the exponent 



sr''' = ---^±y-^ — - — -. (cii) 



which does incorporate correctly the ideal hydrodynamic limit (—4/3) and the Navier-Stokes correction 4/(9ko) but 

is in general higher, the smaller the kq. than the complete IS result S+ = —4/3 + 4/(9ko) — 8/(27ko) + 0{1/kq). 
Therefore, the 'naive' approach overestimates the pressure. In addition, it underestimates the asymptotic pressure 
anisotropy i?^'"^ = 1 — 2/k;o — 8/(3ko) + 0{1/kq), and therefore, overpredicts the magnitude of the shear stress to 
pressure ratio |^|. 



2. Solutions for ultra-relativistic gas with cr2-^2 oc 



A constant cross section implies the existence of some external scale in the problem. For a scale-invariant system, 
however, the only scale available (in thermal and chemical equilibrium) is the temperature, and therefore the cross 
section behaves as ct oc 1/T^. (37), (33) and (29) then give 



(C13) 



i.e., even without the 7r| term, the equations of motion become nonlinear (but are easy to solve numerically). 
For ideal hydrodynamic evolution, p oc t~^^^ and thus, unlike for the case of a constant cross section, 

K{t) = Kof^/^ (C14) 

increases with increasing r. K{t) must grow in general in the viscous hydrodynamic case as well because dissipative 
corrections, namely the ttl/t term in (35), are assumed to be small (or else hydrodynamics is not applicable any 
longer). Consequently, the system gets closer and closer to ideal hydrodynamic behavior as time evolves (as long as 
the expansion is only one-dimensional). For example, the pressure anisotropy approaches unity at late times, for any 
Ko > and initial ttl,o/po, 



The exact Navier-Stokes solution 



Rp{T- 

P™W^(^)'" 



cxd) 



1 . 



1 



4 
3ko 



2/3 



1 



behaves similarly. At late times p oz t "^/^ as in the ideal case, therefore, k{t 
The rate of approach to unity is controlled by the viscosity 

<^(r) = 1 - A (^)'/^ [i + 0{1/kI) + 0((ro/r)2/3) 

Viscosity also increases the pressure relative to the ideal case 



OO) = Ko/iP^f^) - 

_ _2_ru /rn\2/3 
ToTo n 



(C15) 
(C16) 

OO, i.e., iioo = 1- 
(C17) 
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3. Solutions for ultra-relativistic gas with <T2->2 oc t"^^^ 

Near the ideal hydro hmit (i.e., for small viscosities and ttl,o/po), one may substitute the approximate result (C14) 
in the equations of motion (35)- (36) directly. Provided we drop the 7r£ term, these can be converted to a second-order 
linear differential equation, e.g., for p{t): 



11 



40p 2K{t) ( . 4p 



3r 



, 



4po + TTL.O 

3ro 



with initial conditions 

p{to) = Po , p{to) 
The general solution with K{t) from (C14) is[51] 



f-4/3 



p(f) 

7rL(f) = -3f 



.-1/3 ^[T^/^P(r)] 



where 



F±{x) = iFi(±a, 1 ± 2a; -a;) 



(C19) 
(C20) 

(C21) 

(C22) 

(C23) 



are shorthands for confluent hypergeometric functions of the first kind, while C± are matched [52] to the initial 
conditions (C20) 



C± = ±^ boGq:(Ko) - 7rz,,o-FV(Ko)] 



G±{x) 



±2a 



l±2a 



iFi(l ± a, 2 ± 2a, -x) - iFi{±a, 1 ± 2a, -x) 



(C24) 
(C25) 



A very practical approximate formula for the pressure evolution is given by (49), which comes from the asymptotic 
forms (cf. (13.5.1) in [40]) 



iFi(a,6;-a;) 



^^^^ ■ x-" S{a, l + a-b,x) + e^^ (-a:)""'' S{b - a, 1 - a, -x) 



r{b - a) 



Via) 



where 



S{c, d,x) = l + — + c(c+l)ci(rf+l) ^ c{c+l){c + 2)d{d+l){d + 2) ^ 



2!a;2 



3!a;3 



(C26) 



(C27) 



Note that the e ^ term in (C26) is crucial. For large kq, C± are exponentially large, however, the e'^° factors drop 
out [53] in linear combinations relevant for the pressure and shear stress. 

At late times r ^ tq/kq^^ the IS solutions recover ideal hydrodynamics for any initial condition. 



p{t) oc (J^^ ^ , 7rz,(T) oc (J^^ 



To \ 2/3 



for T» 



To 
3/2 



as can be inferred from (C26). The Navier-Stokes solution 

'To\4/3 



P^^(r)=Po(- 



exp 



2 
3ko 



2/3' 



(C28) 



(C29) 



exhibits the same features (as the Reader can easily verify). For the late-time evolution, this scenario gives smaller 
viscous corrections to the pressure and the pressure anisotropy than a oc 1/T^. However, in the large Kq limit we 
recover the same results (C17) and (C18). 
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Analogous derivation gives the exact solutions in the 'naive' Israel-Stewart case: 

p(f ) = C'_ f-2(i+«')/3 f'_{kot^'^) + C'+ f-2(i-')/3 F;(«of'/') , «' = \/f (C30) 



OT 



(C31) 



where 



= ±^ bo G^(«o) - ^0 F4(«o)] (C32) 

= ± a', 1 ± 2a', -x) , G'^{x) = 2x^^F;^{x) - 2(1 ± a')F;^{x) . (C33) 

With the help of (C26) it is straightforward (but somewhat lengthy) to determine the late-time behavior 

^ -r(f)[P(Ko)+^oX(Ko)] (C34) 



Pid. 

where in the 'naive' case 



P— («o) = 1 + — + — + ^ 73 ' ^~M = -;t--^ + -3) • (C36) 



— + — + o(- 

3ko Skq \'«oy ' ^ 2ko 6kq 

Comparing to the 'complete' Israel-Stewart result (C21) (obtained in the small ^ limit) 

we see that for the 'naive' approximation the evolution approaches ideal hydrodynamic p/pideai ~ const behavior 
later (deviation of T from unity is larger), and for near-equilibrium initial conditions (^0 ~ 0) the pressure saturates 
at a higher yalue {P is larger). 

4. Solutions for ultra-relativistic gas with 2 — > 2 cross section and rfa/seq — const 

The last scenario we consider is when the cross section is dynamically adjusted to maintain a constant shear viscosity 
to equilibrium entropy density ratio r]s/seq, such as the conjectured lower bound of l/(47r). From (29), (33), (51), 
and (52), 



1 / ln\f^f(f)]\ , ^ 

^ 1 + \ , C39 



and thus 

Vs Vs,o P{f)f'^ 4-xo 



Seq Seqfi k{f) 4-xo + ln[f4p3(^)j ' 

where 

Seqfi «;o(4 - Xo) 

Therefore, ris/s^q = const requires 



(C40) 



(C41) 



4-xo + ln[T4p3(r)] 
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Within the generic region of validity for viscous hydrodynamics, Itt^I <C p, this scenario also implies a growing 
K(t) ^ T*^^/^ and therefore convergence to the ideal limit at late times. Note that the double ratio (?7s/seg)/ ('7s, o / Seq,o) 
as a function of t/tq depends only on nL,o/Po, i^o, the type of cross section (encoded in K), and xo- 

We now analyze the time evolution of rjg/ s^.q in the three earlier scenarios. Compared to the entropy density, ris/s^q 
contains an additional multiplicative term that comes from the time evolution of the shear viscosity. Assume first, 
for simplicity, that we are very close to the ideal hydro limit, in which case rjs/seq oc t"^/^ /K{t). For a constant 
cross section, this results in a growing rjs/seq oc r^/^; while for the other two cases, a oc r^/^ or a oc we obtain 

"Ha/Seq ~ COnst. 

In reality, there are of course viscous effects. Because 

...^ ^2 C p(f)f^/^ X f^/^ ior a = const 

= \ [p(f)f4/3]3 foraocl/r2 , (C43) 



K{f) 



p{f)f^/^ for a oc r2/3 



the relevant quantity that determines the evolution of rjs/ s^q is pf^/'^. The last term in (C40) is only a logarithm. 
Therefore, the first term, (C43), dominates the behavior. Typically, ttl < and thus dissipation generates an 
increasing pf^/^. The increase in r]s/seq is then fastest for the constant cross section case. The other two cases, 
CT cx and cr cx r^/^, are not equivalent when there is dissipation because for the latter the prefactor (C43) is only 
linear in p(f)f"'/^ and, therefore, rjs/ s^q grows much slower. 



5. Comparison of the various cross section scenarios 

After exploring the general behavior, we compare numerical solutions for the four scenarios. Unless stated otherwise, 
for the r]s/seq = const case we start the evolution from chemical equilibrium, i.e., take xo = 0. For the other three 
scenarios, the pressure and shear stress evolution does not depend on xo- For simplicity, wc start the evolution from 
7rL(7o) = 0, and consider two extremes = 1, i.e., equal expansion and scattering timescales, and Kq = 6.67, i.e., 
6.67 times slower expansion than the timescale for scattering. On all figures, the dotted curves correspond to the 
approximation when the tt\ term in (36) is ignored. 

Figure 8 shows the evolution of the pressure relative to the ideal hydrodynamic p ~ t~^/^ result (for a comparison of 
the same observable between hydrodynamics and transport see Figure 4). Dissipation increases the pressure because 
it reduces the pdV work. The effect is largest for the a = const scenario, while smallest for rjs/seq = const and 
a oc T^/^, which two give basically the same result. For Kq = 1 the fourth scenario a oc is in between these limits 
but by Ko = 6.67 it becomes equivalent to cr a r^/^. Dropping 7r| terms in (36) (thin dotted lines) is a fair 10 — 15% 
approximation for a = const and oc 1/T^ at Kq = 1, which improves to an essentially exact one by Kq = 6.67. For 
the other two scenarios, rjs/sgq = const and a oc t^^^, the nonlinear term can be safely ignored already for Kq = 1. 
Note that for Kq = 6.67, dissipative corrections to the pressure are still very modest 10 — 15% at late t/tq 10 — 20 
in all four cases studied. 

Now we turn to the evolution of the viscous stress ttl shown in Fig. 9. All four scenarios give very similar results 
for the early t/tq < 1.5 — 2 growth in magnitude but they differ in late-time relaxation. As inferred from the 
pressure evolution already, rjs/ s^q = const and cr oc t^^"^ are largely identical and relax quickly toward the ideal limit. 
a — const is the one that stays furthest away from equilibrium. For low Kq = 1, the cr oc 1/T^ case lies in between 
but by Kq = 6.67 it becomes identical to r/s/seq = const and cr oc r^^^. The 7r£ term in the equation of motion 
affects the pressure and the viscous stress similarly, and can be ignored for Kq = 6.67 in all cases - for cr a r^/^ and 
Vs/seq = const cvcu at Kq = 1. 

The same; observations carry over to the pressure anisotropy Rp = pl/pt shown in Figure 10. We plot this quantity 
because it is the same one shown in Figure 1 for the hydro-transport comparison in Sec. V (but note the logarithmic 
time axis there). These results further confirm that cr oc r^^^ is a very good approximation to rjs/seq = const already 
for Kq = 1. 

Figure 11 shows entropy production dS/drj as a function of proper time for the four scenarios, with local thermal 
(^0 = 0) and chemical (xo = 0) equilibrium initial conditions. Due to scalings, only entropy relative to the initial one 
plays a role 



ids Id,) 1 



{dSQ/dm) 4 - Xo 



In(f^p^) 



16 



(C44) 



For Kq = 1, a constant cross section generates about 35% extra entropy by late t ^ 15 — 20ro. With cr oc l/T^, 
the increase is only ~ 30%, while a oc r^/"^ and -qs/seq = const give the smallest increase of about 25%. For a 
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FIG. 8: Pressure evolution from viscous hydrodynamics relative to the ideal hydrodynamic p = Po{to/t) result in 0+lD 
Bjorken geometry for an ultrarelativistic gas with 2 — > 2 interactions. Four scenarios are compared for iC'o = 1 (left) and 5 
(right): a — const (dash-dot-dot), a oc (long dash), a cx r^''"^ (short dash), and rjs/seq = const (solid). Approximate 

results with dropping 7r|, terms in the equation of motion are also shown (thin dotted lines). 




FIG. 9: Same as Fig. 8 but for the longitudinal viscous shear ttl normalized by the initial pressure. 



larger Kq ^ 7, the system is much closer to ideal hydrodynamics and therefore entropy generation is slower - about 
10% for a — const, while only 5% for the other three cases. Note that these results also depend on xo^ but almost 
entirely through the explicit 1/(4 — xo) factor in (C44). Therefore, results for arbitrary Xo 7^ can be obtained via 
straightforward rescaling. In the rja/ Seq = const case the shear stress and pressure evolution also depend on xo but 
only very weakly as we show later below (cf. Figure 13). 

Figure 12 shows the evolution of the shear viscosity to equilibrium entropy density ratio rjs/ S(.q, normalized by the 
initial value of the ratio. The entropy is calculated for a system starting from chemical equilibrium (xo = 0). The 
rough expectations that rjs/s^q ^ t^I"^ for a constant cross section, while rjs/seq ^ const for both a cx 1/T^ and 
a cx T^/'^, hold within a factor of three already for Kq = 1 and up to r = 20to (note that the r^/'^ growth in the 
a = const case has been scaled out in the plots). Relative to this "zeroth order" behavior, for all three scenarios, 
"qs/seq grows with time, reinforcing the general results in Sec. C4. The relative growth decreases with increasing K^). 
The i^o dependence is strongest for the constant cross section scenario: the factor of three gain by r = 20to for Kq = 1 
is tamed to an about 25% increase for ~ 7. For the other two scenarios, a cx l/T^ and a cx r^/'', the ratio stays 
nearly constant much more robustly. As expected (cf. end of Sec. C4), of all cases studied a cx t^I'^ approximates 
Vs/seq = const the best, with only '^10% deviation accumulated by late t = 20ro even for a small Kq = 1. 

Finally, in Figure 13 we show that the results for rjs/seq = const depend only weakly on the initial density, i.e., xo- 
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FIG. 10: Same as Fig. 8 but for the pressure anisotropy evolution 
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FIG. 11: Same as Fig. 8 but for the normalized entropy per unit rapidity [dS / dr]) / [dSo / drj) . 
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FIG. 12: Same as Fig. 8 but for the shear viscosity to equilibrium entropy density ratio rjs/seq. The results for a = const are 
divided by (t/tq)'^^^, otherwise they would quickly grow off the plot. 
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FIG. 13: Initial density dependence of Israel-Stewart viscous hydrodynamic solutions for an ultrarelativistic gas expanding 
longitudinally in 0+lD Bjorken geometry with 2^2 interactions that maintain r]s/ s^q = const. To amplify density effects, the 
initial expansion timescale to mean free path ratio is chosen to be low, Kq ~ 1. Three different initial densities are considered: 
chemical equilibrium n = n^q (solid), oversaturation at n = 4neq (long dash), and undersaturation at n = Ueq/A (short dash). 
Left: time evolution of the pressure anisotropy Rp = pl/pt- Right: time evolution of the produced entropy per unit rapidity, 
normalized to the initial entropy per unit rapidity. The produced entropy is scaled by (4 — Xo)/4 to eliminate trivial density 
effects that do not come from the shear stress and pressure evolution (see text). Approximate results with dropping 7r|, terms 
in the equation of motion are also shown (thin dotted lines). 



Density dependence in shear stress and pressure evolution arises in this case because the cross section is a function 
of the initial density (see (C42)). The dependence is weaker, the closer the system is to ideal hydrodynamics because 
in that case p oc r"'*/"^ and xo drops out from K{t). But even for a pessimistic Kq = 1, the pressure anisotropy (left 
plot), varies less than 10% as we change the density by a factor of 4 around chemical equilibrium density {xo = ± ln4). 
In fact a decrease in the density has a much weaker effect than an increase. The right plot shows the effect of the 
same initial density variation on entropy dS/drj production normalized to the initial entropy. Most of the density 
dependence in the entropy change comes from the trivial 1/(4 — xo) prefactor in (C44) which is there in any cross 
section scenario even if the shear stress and pressure evolution are independent of the density. To highlight dynamical 
density effects, we therefore plot, again for a pessimistic Kq = 1, the normalized change in entropy 

4 - xo AidS/drj) _i-Xo f (dS/drj) _ \ ^^^^^ 



4 {dSo/dv) 4 \{dSo/d7^) 

(the scaling factor is chosen such that it has no effect for chemical equilibrium initial conditions xa = 0). The results 
show practically no density dependence, apart from few-percent changes, even for such a low Kq. 

APPENDIX D: USEFUL RELATIONS FROM COVARIANT TRANSPORT 
1. Particle number and transverse energy 

The particle number and transverse energy distributions for particles crossing a 3D hypersurface cr{x) = const are 
given by 

dN = dydp'^j, J p^da^(a;)/(x,p) (Dl) 
dEr = dydp^ J pt'daf,{x)mTf{x,p) (D2) 



where m-p = yjp^ + , Pt = yPx ^ Py transverse momentum, and da^^{x) is the normal to the hypersurface 

at space time coordinate x. For our boost-invariant scenario it is natural to follow quantities per unit coordinate 
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rapidity as a function of the proper time r. For r = const hypersurfaces p^da^ = rriT rchw (PxTdrj, and in our 0+lD 
case, / only depends on shw, p±, and r, where co = y — rj. Thus, 



dNjr) 

dr] 
dExir) 
dr] 



= rArld^Prd.mrc^.fir,sl^,Pr) 

= tAt I d^ pt dio nix chijj f{T,sh.ui,pT) ■ 



(D3) 

(D4) 



At is the transverse area of the system. With the local thermal equilibrium distribution for ultrarelativistic particles 



87rT3 



and the quadratic form (62), straightforward integration gives 

dN_ 
dr] 

dErir) 



n T At = const 



dr] 



SnTdN / 5^ 
~~dr] V ~ 16 



(D5) 

(D6) 
(D7) 



Clearly, dissipation slows the decrease of transverse energy (for typical ttl < 0), and 2 — > 2 interactions of course 
course conserve particle number. 

Note that dET/dr]/{TAT) is almost identical to the transverse pressure (DIO), but has an extra chw factor in the 
integrand. 



2. Early pressure evolution 

Here we evaluate the early transverse and longitudinal pressure evolution from the transport for a local equilibrium 
initial condition. The results hold for any interaction, not only 2^2. 

In local equilibrium the collision term vanishes, thus in the vicinity of r = tq the evolution is governed by free 
streaming. In our 0+lD case, free streaming 



chujdr — ——da. 

T 



/(shw,p±,T) = 



{uj = y-r]) 



implies 



/(sha;,p_L,r) = /(Tsha;/ro,p_L, tq) . 



(D8) 



(D9) 



Substituting a local thermal initial distribution for ultrarelativistic particles (D5), the definition of the energy- 
momentum tensor 



T^-'iV = 0, r) = / ^p^p- f=[ d^pi, dyp^p- f{shy,p^, r] 
J Po J 

gives the transverse pressure 

Pt(t) = T^^(j7 = 0,t) = Af J dp±p±d^dy {p± cos (1))'^ exp 
Here a = t/tq. Change of variables to q = ashy leads to 



(DIO) 



— Vl + a^sh y 



dy 



3Ton 

J (1 + a2sh2y)2 



(Dll) 



Pt(t) 



3Ton 



oc 

/ 



dq 



= Ton 



4(a2 - 1)3/2 



acos- 



Analogous calculation gives for the longitudinal pressure 



PL{T) = T''{r] = 0,T) = 3Tor 



dy sh y 



(1 + a^sh^yy 



Ton 



2(a2 - 1) 



acos- 1 



(D12) 



(D13) 



29 



Expanding near a = 1 



Pt{t) 



Pl{t) 



Ton 1 



Ton 1 



4(r - To) _ 

5to 
12(r - To) 

5ro 



+ 0((t-to)^) 



+ O{{r-ror) 



(D14) 



(D15) 



and thus (65) follows. 
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